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Abstract. We calculate the discrete moments of the characteristic polyno- 
mial of a random unitary matrix, evaluated a small distance away from an 
eigenangle. Such results allow us to make conjectures about similar moments 
for the Riemann zeta function, and provide a uniform approach to understand- 
ing moments of the zeta function and its derivative. 



1. Introduction 

Let £(s) be the Riemann zeta function, and denote its non-trivial zeros by i + i7n f 
with < 71 < 72 < . . . . (For simplicity, we will assume the Riemann Hypothesis, 
which says that j n el). It is known (see, for example, Titchmarsh's book Q for 
details) that if N(T) is the number of zeros with < 7„ < T then 

iY(T)= £ iog 2^ +o(iogT) - 



Moments of the zeta function, 



T 

2k 



it. 



have long been of interest to number theorists, with it being widely believed that 
1 



T 

with 



T \(a + it)\ 2k dt ~ f k a(k)(\ogT) k 



o 



p prime 



and fk being an integer (when k is integer) divided by (fc 2 )!, whose value was 
unknown apart from in a few cases. The known values of /& are fo = 1 (trivial), 
fi = 1 (Hardy and Littlewood §]), and f 2 = ± (Ingham g). The value fs = § 
has been conjectured by Conrey and Ghosh Q and = 2 ^" 24 is a conjecture of 
Conrey and Gonek || . We should mention that a(k) given in (L.l) can be calculated 
for certain k: a(0) = a(l) = 1 and a(— 1) = a(2) = 

In H Keating and Snaith argued that one can create a probabilistic model for the 
zeta function around height T using the characteristic polynomial of an N x N 
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unitary matrix chosen according to Haar measure, when 

N = l °^- 

Setting 

Z v {9) := dct{I~e- i6 U) 

N 

= n( 1 _ e K^)) ) 

n=l 

and defining Mjv(2A;) := E^r {|Z;/(0)| 2 ' C } where Ejv denotes expectation with re- 
spect to Haar measure, they found that 

G 2 (fc + 1) fc * 
~ G(2fc + 1) 

as N — > oo for fixed fc subject to 9ie(fc) > —1/2, where G(-) is the Barnes G- 
function. 

By comparing with the known (and previously conjectured) values of they were 
led to conjecture that = Q/2fc+i) • 

Conjecture 1. (Keating and Snaith). For fixed k > —1/2, as T — > oo, 
1 r^ix-^* G 2 (fc + 1) r^A r V 2 

where G(-) is the Barnes G-function and a{k) is given in (1.1) 

Following this, Hughes, Keating and O'Connell [Q used the characteristic polyno- 
mial to model the discrete moments of the derivative of the zeta function, 

jk{T):= wD £ icti+^r- 

V ' 0<7„<T 

They calculated that for fixed k subject to 9te(fc) > —3/2, 



as N — > 00, and they used this to conjecture the asymptotic form Jfc(T) should 
take for large T. 

Conjecture 2. (Hughes, Keating and O'Connell). If all the zeros of the zeta 
function are all simple, then for fixed k > —3/2, as T — > oo, 

1 iw/i • M 2fc G 2 (fc + 2) „./ r\ fc(fc+2) 

^ o E <T lc'(l+^)l -G(2TT3y a ^( log ^ 



where a(k) is given by (1.1), and G(-) is f/ie Barnes G-function. 
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Again this conjecture is found to agree with all previously known results; when 
k = — 1 (a conjecture of Gonek ||), when k = (trivial), and when k = 1 (a theo- 
rem of Gonek H under RH). Also, extending a theorem due to Conrey, Ghosh and 
Gonek ]To| (recorded in this paper as theorem^) beyond its range of (proven) appli- 
cability, conjecture [Tl] states that J2(T) ~ 144 1 07r ^ (log , which agrees perfectly 
with conjecture |^. 

It is striking that conjecture |l| and conjecture ^ have very similar form. The purpose 
of this paper is to unify them as special cases of one result. 



In the next section we will show 



2/r 



G 2 (fc + 1) 



F k (y)N k 



G(2k+iy 

where F k (y) is a certain function, independent of N, given in theorem | 
We will then use this to conjecture that for L = ^- log 



1 



T 



]VCf) o< E t IC (I + i (7. + a/L)) \ 2k ~ ^±^ a (k)F k (2na) (log - 

and will show, in section [|, that this conjecture contains conjectures [j] and |^ as 
special cases (a — > oo and a — > respectively). The conjecture is found to agree 
with a known result of Gonek ||] , and the extension of the theorem due to Conrey. 
Ghosh and Gonek cited above. These comparisons are discussed in section 3T 



2. The random matrix calculation 



Theorem 3. For fixed k with f Kt(k) > —1/2, and for x < AN with A < ir constant, 



2x 
TV 



2k 



G 2 {k + l) 



F k (2x)N K +0[N 



rk'-X 



J G(2fc + 1)' 
where 

(2.1) F k {2x) = x 2 ]k {x) 2 + x 2 j k -i[xf - 2kxj k (x)j k -i(x) 

where j n {x) are the spherical Bessel junctions of the first kind. 



Proof. First note that 



\Z v (fi N +0) 



21, 



1 



-i/3 



2 k 



N-l 







1 -0n-I3) 



n=l 



2 k 



The average of this over all N x TV unitary matrices with Haar measure can be 
written (JO, O) as an TV-fold integral 




N-l 



EH 



n = l 



2k 



Mi... ae 



N 
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Putting all the j — N terms from the first product into the second, we see that 

(2.2) E N [\Zu (8 N + [3)\ 2k } = 1 |2sin (\(3) \ 2k Ejv-i | |Z^(0)| 2 \Z v {(3)\ 2k ) , 

where is the characteristic polynomial of an [N — 1) x [N — 1) unitary matrix. 
By rotation invariance of Haar measure, 

Ejv-1 {|^(0)| 2 |%(/3)| 2fc } = Ejv-i { |^(0)| 2?C | V^T} 
This is calculated in Theorem^ below, where it is shown that for 9 c te(fc) > —1/2, 

E iV {|^(o)| 2fc |^ t/ (y/iV)| 2 } 

G 2 (k + 1) ~ k(k-l+p)\(k + P )\ 2 (fe+1)2 

G(2/c + 1) ^ p!(2fc + p)\(2k + 1 + 2p)\ 1 J tf 



and substituting this into (2.2) (where we put (3 = y/N) we see that 



E 



Zu [On + 



21, 



G 2 {k + l) 
G(2k + 1) 



F k (y)N k 



where 
(2.3) 



Flc(y) - k S^ (k-l+p)\(k+p)\ )p2k+2p 
tk( - y) - k 2^p\(2k+p)\{2k + \ + 2p)\ { l) V 



The spherical Bessel functions of the first kind are defined as 



jn{z) = J —Jn+l/2(z) 

^ m!(2n + 2m+ l)!2 2ml j 
where Ji,(z) is the z/-th order Bessel function of the first kind. Hence, 

F k (2x) = x 2 j k (x) 2 + x 2 j k ^ 1 (x) 2 - 2kxj k (x)j k ^ 1 (x) 
which can be seen by comparing the Taylor expansions. 

The above large- N asymptotics are for x = o(N). This can be extended to \x\ < AN 
for A < 7r an arbitrary constant as follows: Let (3 be a fixed constant subject to 
< (3 < 2tt. By (|J) and the results of Basor |||, 



E 



N {\Z U (6 1 +(3)\ 2k ] = 1 ^sinQ/?)! 2 ^^ {|^(0)| 2fe \Z D ((3)\ 2 } 
G 2 (k + 1) 



G{2k + 1)' 

If one lets x = N/3/2 then the large- a; asymptotics of Bessel functions (see, for 
example, chapter 9 of [[l4|) implies that 

1 



F fc (iV/3) = 1 + 



and so theorem^ gives the correct first order term as N — > oo, for \x\ < AN with 
A < 7r a constant. □ 
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j n (x) = (-1)"*" 



Remark. When n is an integer, 

1 d \ " sin x 
v x dx ) x 

which leads to a neat evaluation of Fk(2x) for integer k, the first few being: 

x 2 — sin 2 (cc) 



Fi(2x) = 

, . „ ,„ , x 4 -3x 2 +3a;sin(2;r) + (2a; 2 -3) sin 2 (a;) 
2.4) F 2 (2x) = — i '- ^ 

. , x 6 -3a; 4 -45a; 2 + (-12a; 3 + 45a;)sin(2a;) + (-3a; 4 + 72a; 2 -45)sin 2 (x) 

F ^ 2x) = ^ 

Theorem 4. 

EAr{|Z(0)| 2fc |Z(/3)| 2 } = M N (2k)Nl{N + 2k)lx 

x " (2sin(/3/2)) 2 » (JV + A:-m)!(fc + m + n)! cW2m _ Ar+n) 
^ n\(2k + n)l ^ (N-n-m)\m\ 

n=0 v ' m=0 v ; 

where M n (2k) is given in (1.2). If & — ► as N ^ oo, then for $He(fc) > —1/2, 

e*(|z(0)| |z(-)| } = G(2fc - 1} ^ p!(2fc + p)!(2fc + J + 2p)! 



Proof. If k in an integer, then 



E 1 »mjl-e""| 2 *|l-e" , -e-"'| 2 J 



_7r l<j<m<JV 

JJ(1 - e ie ") fe (l - e - w ») fc (l ~ e ie "e~ i/3 )(l - e - ie "e i/3 ) d0„ 



JV 

X 

n=l 

which equals, after some simple manipulation of the terms 



(_\\N{N-l)/2+kN+N e -iN0 

(2tt) n N\ X 



J... J Y[ (e ie > -e w '») 2 JJ(e ie »)- JV - fc (l-e w -) a *(e ifl »-e i ^) a (W B 



_7r l<j<m<N 

Therefore, 

"TV! 



E^IZ^O)! 2 *^/?)! 2 } = £_ ( _i)*C"-i)/a+**+" 



CT i J] ft; - t m ) 2 n - t«) 2fe (*n - e i/3 ) 2 1 

[ l<j<m<N n=l " J 
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where CT{-} denotes the constant term in the Laurent expansion in the variables 
ti, . . . ,tjv- The constant term equals (by lemma 1 of fl5|| ) 

r r 1 N 

V ^ •* ^° l<j<m<N n=l 

Kaneko has evaluated this integral (which is a generalization of Selberg's inte- 
gral) as 

T(l+j)T(j + y-N- k)T(j + 2fc + 1) 
11 T (j+y + k + l) 

(-N) m+n (y + k + l) m+n e 2 "^(l - e^) 2 " 
J^=o (y-N-k) m (2k + l) n mini 

where (a)„ = a(a + 1) . . . (a + n — 1) = T(a + n)/T(a). Since we have assumed that 
k is an integer, 

f_l\N+k-j 

lim t/r(t/ + 7 — N — k) = — ; r 

y^o J u J ' T(N + k-j + 1) 

and so we have 



E^{|^(0)| 2fc |Z t7 (/3)| 2 } = ^-n 



li rO' + * + i)r(JV + A-j + i) 

(-JV) m+n (fc + l) m+n e 2 ""^(l - e^) 2 " 
^ n (-N-k) m (2k+l) n mini 

m,n— 



Expanding everything out in terms of the gamma function: 

e n {\ Zu m 2k \zum 2 } = M N (2k) T(N + 1 + 2k)m 



r(N + i + fc) 2 

A (2sin(/3/2)) 2 " y? r(JV + fc + 1 - m)r(fc + n + 1 + to) e^^H 
^ n!T(2fc + n + l) ^ T(iV + 1 - n - m)\ ml 

n— m— 



where M^(2k) is defined in ([L4). Observe that the inner summand is invariant as 
to — ► N — n — m, and so the inner sum is in fact a sum of cosines (and thus the 
series expansion in contains only even powers of [3). Furthermore observe that 
both sides of the equation are analytic functions of k (for 9le(fc) > 0), both sides 
can be easily bounded bounded by 0(2 2N ^ t ^ k ') (for large k, with N and (3 fixed), 
and the two sides are equal at the positive integers. Thus by Carlson's theorem (see 
§17 of p2[ , for example), the restriction on k being an integer is no longer required, 
and the left-hand side equals the right-hand side for all complex k. 

Now, using the fact that (— a) n = (— l) n T(a + 1)/T(a + 1 — n), we have 

y? T(N + k + 1 - m)T(k + n + 1 + to) e W 2 m-N+n) 
^ T(N+l — n-m) ml 

m=0 y ' 

= r(iV + fc + l)r(fc + n + 1) ^ ( -N + n, k + n + 1 _ p2ip \ pi{n - N) p 



T(N -n+1) V -N-k 
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where 

'a,b \ _ {a) m (b) m z r< 



c J ^ (c) m m\ 

/ m—Q v ' 

is the Gauss hypergeometric function. (If a, c are negative integers with —a < — c 
then it is a polynomial of degree —a). 

Applying the quadratic hypergeometric transformation 15.3.26 of (m) we get 

(2.5) 2 F 1 (- N + n ; k+ i n + 1 -e^)e^ N ^ 
\ —N — k J 

For m a positive integer, 

/ —771,6 \ r(l — C — 77l)r(l — C + 6) / — 771, 6 

21 v c ; 7 = r(l-c)r(l-c-m + 6) 2 1 v l-c-m + 6 ;1_Z 
and so we see that the right-hand side of ( [2.5|) equals 

r( ijv + |n + fc + i)r(|iv + !« + * + §) n> 



r(iv + /c + i)r(/c + n + f) 



-(2cos/3) JV - n x 



/ -IiV + In,-iiV+i77+i 

X 2-^1 , a ,1 - 



77 + k + § cos 2 (/3) J 

Therefore, we have proven 

N l± 2 (N-n)\ 

T{N,k,m,n)(-l)"\2sm(f3/2)) 2n (smf3) 2m (cos(3) N - n - 2m 

n—0 m— 

where 

T(N,k,m,n) = M N (2k) 



Nl(N + 2fc)!(iiV + in + fc)!(iJV + \n + k + i)! 



(A^-77)!(iV + fc)! 2 
2 JV ""(-iiV + |n) m (-iiV + in + i) m (fc + 77)! 



m\n\(2k + n)!(n + fc + | + 777)! 
Observe that for fixed k,m,n with 9te(fc) > —1/2, 

Tf/V ft 777 n) = ^(fc + l) (fc + 77)!(fc + 77 + 777)! 

V ' ' ' ' G(2k + l)m\n\(2k + ny.{2k + 2n + 2m+iy. 
and so 

{|Z(0)| 2fe |Z(t,/A0| 2 } ~ g^ll /Y^ 1 ) 2 



OO OO 



(fc + n)!(fc + n + m)! „, 2 m+2n 



J m!n! (2/c + n)!(2fc + 2n + 2m + 1)! 

n—0 m=Q 



s 
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in the sense that for each fixed integer h, the coefficient of y h on the Taylor expan- 
sion of the left-hand side converges to that of the right-hand side as N — > oo. 

Finally we show that 



oo oo 

EEC- 1 ) 

n— m= 



(fc + n)\(k + n + m)! 



2m+2n 



m!n!(2fc + n)!(2fc + 2n + 2m 

OO 

=E 



fc(fc - 1 +p)!(/c 



^ Q p!(2fc+p)!(2fc + l + 2p) 



■(-i)V 



This can be proven by comparing the coefficients of y 2p . That is, we wish to show 
for all integer p > 0, 
v 



E 



(-l)P~ n (k + n)\{k + p)! 



(-l)Pfc(fc- l+p)!(fe + p)! 



^ (p - n)!n!(2fc + n)\(2k + 2p + 1)! 



p!(2fc+p)!(2fc 
This is equivalent to showing 

A {-l) n {k + n)\p\{2k+p)\ 



l + 2p)\ 



(2.6) 



k{p - n)\n\(2k + n)\(k 



= 1 



l+p)\ 

for all integer p > 0, and this we shall do by creating the Wilf-Zeilberger pair jlj 
Denote the summand in (2.6) by F(p,n), and observe that 

(2.7) F(p + 1, n) - F(p, n) = G(p, n + 1) - G(p, n) 

where 



G(p,n) 



(2fc + n)ri 



-F(p, n) 



(n - p - l)(fc + p) 

(the (2/c + n)n/(n — p — l) (k + p) being calculated by Zeilberger's algorithm). 
Summing both sides of (2.7) over n, we see that the right-hand side telescopes to 
zero, which shows that the left-hand side of (2.6) must be a constant, independent 
of p. Putting p = 0, direct calculation shows that constant is 1. □ 



3. Conjecture about the zeta function 
Conjecture 5. For fixed k subject to D\z(k) > —1/2, 

as T — > co, uniformly in a for \a\ < L, where L = j- log w- is the density of zeros 
of height T. G(-) is the Barnes G-function, a(k) is given by (1.1), and Fk(2ira) is 
given in theorem 

If this conjecture is true, then we are able to prove conjecture || and a variant of 
the Keating-Snaith conjecture (conjecture |l]). 

Corollary 5.1. If conjecture^ is true, then 

1 ,„n . M2fe G\k + 2) / T\ k{k+2) 



N(T) o<f<T G(2k + 3) "\ *2n 
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Proof. By the definition of differentiation, 

1 2k 



| C '(|+i7 n )r = i 2fc lim 
From ( |2.3|) we have 



K(| + i(7« + f))| 



,. F(2na) ,„ s2k k\k\ 



a^O a 2k v ; (2fc)!(2fc + l)! 
so applying conjecture |5| and using uniformity to swap the a — > and N — > oo 



limits, we have 



1 V- w/i ■ M2fc G 2 (fc + 1) 2fc (2^) 2fe fc!£;! / T 



0<7„<T 



G 2 (k + 2) ( T ^ k[k+2) 



a(k) log — 



G(2k + 3) w V 27r , 
as required. □ 

Corollary 5.2. From conjecture |^ follows that for (3 > /reed 
<«> iv^ o< E T MI + '(7„ + «)| 2k ~||^a W (lo g ^ ^ 

Proof. The asymptotics as z — ► oo of the spherical Bessel function (see chapters 9 
and 10 of 111) are 



jn (z) 



I(_l)(«+i)/2 C0SZ if n is odd 
-(— l) n / 2 sinz if n is even 

and putting this into ( |2.l| ) we have 

lim F fc (27ra) = 1 

a — >oo 

Setting a = L[3 in conjecture || completes the proof. □ 

Remark. Note that corollary |5.2| can be thought of as a variant of conjecture |[ 
This is because one expects the mean of |C(l/2 + i£)| 2fe to be independent of the 
average taken, and the Keating-Snaith conjecture is a result about the continuous 
mean, whereas corollary [T^ is a result about a discrete mean. To see this, recall 
that the zeros get denser higher up the critical line, and so if (3 > is fixed and j n is 
random, one might expect £ (h + i(7„ + /?)) to be random (whereas, if (3 was small, 
it would be highly influenced by the fact that £(l/2 + i7„) = 0). The left-hand 
side of (3.1) averages this, and thus acts like a discrete mean of |C(l/2 + ii)| 2fc . 

3.1. Comparison with the zeta function. Gonek || showed that if the Riemann 
Hypothesis is true then 

1 iwi ■( it\\\ 2 I i /sin(7ra) x2 



|C(i + i(7» + a/£))| 2 ~fl-(- 




N(T) ^ '^ 2 w " • 1 "\ \ \ na 

0<7„<T \ v 

uniformly in a for \a\ < L/2, which is in perfect agreement with conjecture ^Jwhen 
k = l. 
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There is no proof of the conjecture for k — 2 (unlike conjecture ^ which is proven 
for k = 1 and 2). But there are theorems along the lines of conjecture || for k = 2: 

Theorem 6. (Conrey, Ghosh and Gonek pj.) Assume GRH and let A(s) = 
J2n<x n ~ S where x — for some r\ £ (0, |). TTien, 

7~T > \CA (h + 1 {in + a/L))\ ~^7> T7 x 

N(T) ^ |S 12 Wn ' ;;| 7i- 2 ^ (2j + 5) 

V ; 0<7„<T j=0 w y 



x ^ + |(2j + 5)^-|^^+ 8 + ^+ 7 + ^(l-^+ B ) (log J)' 
uniformly for bounded a. 

(We have slightly changed notation from [JlO|, to be consistent with our definition 

Putting 77 = 1 in the above (which, as it stands, is not allowed under the conditions 
of the theorem) then A{\ + it) = C(| + it) + 0(i~ 1/2 ), and we have 

E \C 2 (h+Hln + a/L))\ 2 

0<7„<T 

(2/+ 6)! (^ + 5 ^( 10 ^ 
Note that 

4 ^ (-1)^(2^)^ 
^Z. (2i + 6) , (2j + 5jj- 

1 . . (2tt 2 q; 2 - 3)sin 2 (7ra) +37rasin(27TQ) + (vra) 4 - 3(7ra) 2 
12 a(2) (^ 

which is what is predicted in conjecture |^. 

That is, from a purely number theoretical calculation involving no random matrix 
theory, we have 

Conjecture 7. Assuming that rj = 1 is permissible in theorem^ then 

^ o< E t lc (\ + i (7« + a/L))\ 4 - ^m^) (log |) 4 



where F<i(flx) is given in (2.4). 



So, following the proof of corollary 5.1, we may deduce 
Corollary 7.1. If conjecture^ is true then 

Jjffi E l^(i +i ^)| 4 ~144^( l0S £ 

0<7„<T v 

Note that this is the same answer that one gets from putting k = 2 into conjecture 
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